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pi • Abstract. Collisions between tightly confined atoms can lead to ionization and hence 

to loss of atoms from the trap. We develop second-order perturbation theory for 
a tensorial perturbation of a spherically symmetric system and the theory is then 
applied to processes mediated by the spin-dipole interaction. Redistribution and loss 
^ , mechanisms are studied for the case of spin-polarized metastable helium atoms and 

' results obtained for the five lowest s states in the trap and trapping frequencies ranging 

• rH ! from 1 kHz to 10 MHz. 
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1. Introduction 

There is significant interest in the study and control of quantum processes involving 
trapped ultracold neutral atoms where the trapping environments are so tight that 
the effect of the trapping fields upon the colliding atoms cannot be ignored or 
approximated as constant background fields. Trapping in three-dimensional optical 
lattices, with typical trapping frequencies of 10^ to 10^ Hz, forms the basis of such 
studies as quantum phase transitions of ®^Rb atoms 1|, storage of metastable argon 
atoms |2], implementation of quantum logic gates and formation of highly entangled 
quantum states [HI Ej. Theoretical investigations have focussed mainly on tightly 
confined alkali systems and have been based either upon direct numerical integration 
of the radial Schrodinger equation for the relative motion of the colliding atoms using 
the best available full interatomic potentials ^ or upon a regularized delta-function 
pseudopotential and an energy-dependent effective scattering length [7] . 

An understanding of collision processes in trapped ultracold metastable 2^S helium 
(denoted by He*) is necessary to obtain Bose- Einstein condensation of this species |H1IS] 
and to investigate these novel excited-state condensates |TIH ITT] . Although current 
experiments on He* only use trapping frequencies of the order of 10^ to 10^ Hz, it is of 
interest to investigate the effects of much tighter trapping on the allowed quantized trap 
states, as a possible tool to manipulate the confined atoms, and to enhance trap loss 
through ionization processes at small interatomic separations as a means of studying 
these processes. 

We have recently analyzed a system of two colliding ultracold atoms under strong 
harmonic confinement in a spherically symmetric trap from the viewpoints of quantum 
defect theory and of elastic scattering in the interatomic potential. We have developed 
methods for determining the energies of the quantized states produced by the presence 
of the trap and the theory was applied to collisions between spin-polarized He* atoms, 
see Peach et al \\2\ IT^ . The energies were determined for a wide range of trapping 
frequencies for s- and d- wave collisions using two totally independent methods to 
integrate the radial Schodinger equation. Excellent agreement was obtained between 
the two methods, one based on the use of quantum defect theory and the second on the 
use of a discrete variable representation. 

These calculations ignored loss processes, but inelastic collisions may cause 
transitions to states from which there is a high probability of Penning and associative 
ionization. A study of such loss processes is the subject of this paper which is organized 
as follows. In section 2, the theory of collisions in an isotropic trap is briefly reviewed, 
and in section 3, second-order perturbation theory is introduced for a general form of 
the perturbation and for trap states of any angular momentum. In section 4, the theory 
is applied to perturbation by the spin-dipole interaction and in section 5 the numerical 
methods are described. Finally in section 6 the theory is applied to the case of spin- 
polarized He* atoms. It is found that only the s states are significantly perturbed and 
shifts and lifetimes are presented for the five lowest s states in the trap for trapping 
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frequencies ranging from 1 kHz to 10 MHz. 



2. Collisions between two atoms in an isotropic harmonic trap 



Consider two atoms with masses Mi and M2, spin quantum numbers 5*1 and 5*2 and 
position vectors ri and r2 relative to the centre of the trap. The interatomic separation 
is given by r = |r| = |ri — r2| and, for the case of an atom- atom potential that is only 
a function of r combined with a potential for an isotropic harmonic trap of angular 
frequency cu, the Hamiltonian is separable into two parts iJcm and Hq describing the 
centre-of-mass and relative motions of the two atoms. Here we will use and extend the 
notation developed in ^Hj- If the total spin quantum number is S and the adiabatic 
potential for the molecular state ^s'+i^ jg denoted by Viys{f)i the equation for the relative 
motion is 



Ho\i^{r)) = E\ij{r)) 
where Hq is defined by 

_o 1 



2M 



trap 



(1) 



(2) 



In ((H) and (j2I) is the energy eigenvalue, the reduced mass M = M1M2/ {Mi + M2) and 
^trap = 1 or according to whether the harmonic potential is turned on or off. If the 
angular momentum quantum number for the relative motion is Z, the eigenvector \ip{r)) 
is given by 



iV'(r)) = -FuiKs{r)\lm)\SiS2SMs) 
r 



(3) 



where the magnetic quantum numbers m and Ms refer to projections of the angular 
and spin momenta onto the molecular axis. The radial function FkiAs{r) satisfies the 
equation 

^ d^ _ /(/ + 1) . 2M, 



■VAs{r) + e 



(4) 



where we have introduced the quantities 



2ME 



h 



(5) 



' ^ Mlu 

It has also been shown in ^Hj that for the discrete states with E > produced by the 
presence of the trap, hereafter referred to as trap states, it is natural to introduce an 
effective quantum number (or scaled energy) n*, where 



n 



E 
2Kuj 



, I 3 



0,1,2, 



(6) 



and fi' is a quantum defect that varies slowly as n'^ increases. By introducing the 
dimensionless variable p = r/C,, equation can also be written in the form 

J 9 9 P ^trap (- 

ap^ p^ fiuj 



An* 



FkiAsip) = 0. 



(7) 
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3. Perturbation theory 

The Hamihonian for the perturbed system is given by 

H = Ho + Hp, (8) 

where Hp is the perturbing potential and it is assumed that the eigenstates of Hq are 
known, i.e. 

Ho\j) = E,\j); J =0,1,2,.... (9) 
Then the change in energy of the initial state i, correct to second order, is given by 

AE = AEi + AE2 , (10) 

where 

AE, = {t\Hp\t) (11) 

and 

The calculation of AEi is straightforward and AE2 can be evaluated as follows. Using 
the method of Dalgarno and Lewis |^ , we introduce the operator F which satisfies the 
inhomogeneous equation 

FH^-H^f\ \{) = Hp\i), (13) 

so that p2|l becomes 

AE2 = ml - {i\Hl\i){i\F\i) . (14) 

This result is only useful if we can determine the operator F. We consider the 
perturbation given by 

i/p = T(A)-U(A)\/p(r), (15) 

where T(A) and U(A) are tensor operators of order A and Vp(r) contains the radial 
dependence of Hp. We set \i) = {ipii^)) and make the following expansion 

F\t)=F^FkiAs{r)\lm)\SMs)= ^ lfi'm's'M',ir)\l'm')\S'M's) , (16) 

I'm'S'M'g 

where we have suppressed the spin quantum numbers 5*1 and 5*2 in Q. Then, on using 
dH)-©, (HSl), (Hi) and (Uni), we obtain 

{FH,-H,m = ^ Y: 



2M 



Vm'S'ML 



X 



d^ /'(/'+ 1) 2M ■ 

Tl^trap ^V^'S' + 



fl'm'S'M'\l'm')\S'M's) 



dr2 r2 ^4 t^^^P 

T(A)-U(A) VpFMAs\lm)\SMs) . (17) 



Collisional loss of atoms from harmonic traps 5 

Now if we define GkVK'S'ij') by the relation 

fvm'S'M',{r) = {S'M's\T{X)\SMs)-{l'm'\{J{X)\lm) Gh'A's'W , (18) 
then Gki'A's'i'^) satisfies the inhomogeneous radial equation 



c.f. On introducing the scaled variable p, ()19p becomes 



Gki'A's'ir) = ^V^p(r) F^iAsir) , (19) 



p Atrap ; h 4n 



2\/p(p) 



Gkl'A'S'ip) = — FkiAsip) 



dp2 p2 ^ — P 

= W^(p)i^fc^A5(p), (20) 

c.f. ((Tj). If Atrap = in (jD), then in the outer region where V^sir) is very small, the 
regular and irregular solutions for FkiAs{f) are given by 

FP'{r) ~ (A;r)[cos5i ji{kr) — sin 5; ni{kr)] ; 

F/(r) ^ -A^ (A;r)[sin5i j'K^^) + cos 5; ni{kr)] , (21) 

where ji{kr) and ni{kr) are spherical Bessel functions, see ^SJ; is a normalization 
constant and 6i = 6i{k) is the /-wave phase shift for elastic scattering. The wave function 
FkiAs{f) is matched to a normalized bound-state wave function of the same energy by 
choosing 

dp'' 

dn* 

c.f. (jni), where p' = p'in*) is treated as a continuous function. Also, as r — > oo, it can 
be shown that the contribution to the solution of ()19|) from the particular integral is 
given by 

Gm'A'S' = "^[[Fhr] + ii^/(r)] ^i?(x)| ; x = ^, (23) 
where R{x) is a slowly varying complex function of x satisfying the conditions 

= . (24) 



A^ 



T^ek 1 



(22) 



^ M dR 
R(0) = — i- — 

^ ' 2k ' dx 



x=0 

In what follows we use the results in (I14j) - ((201), average over initial degenerate 
states \lm) and carry out some angular algebra, more details of which are given in the 
Appendix. We introduce reduced matrix elements (j'||X(A)||j), see and then /S.Ei 
in (fTTj) is given by 



A^i = [(2/ + 1) (25 + l)]-5 {S\\T{\) \\S) (/| |U(A) 1 1/) / FuiKs{r)V^{r)FkiKs{r)dr 

Jo 

(25) 

for A = and zero otherwise. On using (fTTjl . (jA.ljl and ()A.2jl . AE2 in (fT^ becomes 
AE, = [{21 + 1)(2A + l){2S+ir' 5^ |(S'||T(A)||S)p |(n|U(A)||/)p 



I'S' 
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GkVK's' {'r)Vp ir)FkiAs{r) dr 

POO POO 

- 6in5s's / FkiAsir)Vpir)FkiAsir)dr / GkiAsir)FkiAsir)dr . (26) 
Jo Jo 

4. The spin-dipole interaction 

The interaction between the electronic-spin magnetic-dipole moments of each atom 
produces the spin-dipole interaction Hamiltonian 

i^sd = Vpir) 1 [3 (Si ■ r)(S2 ■ r) - S, ■ S^] , (27) 

where Si and S2 are the electronic-spin operators for the two atoms and f is a unit 
vector directed along the internuclear axis. The function Vp{r) is defined by 



/3 ^ 2 / A'-e 



W = -^; P = a'(^{E^al (28) 

where a is the fine structure constant, ao is the Bohr radius, (/ic/ /ie) = 1.00115965 is the 
electron magnetic moment and is the Hartree energy (= 1 a.u.). The perturbation 
-f^sd in can easily be identified with Hp in (fT3j) since 

S1S2 = ^ [SS - SiSi - S2S2] , (29) 

where S = Si + S2 is the operator for the total spin. Therefore in ()15j) . A = 2 and 

T(2) = S{2) - 5i(2) - 52(2) ; U(2) = i(3ff - I) , (30) 

where I is the unit dyadic. In (jHOl), <S{2), iSi(2), 52(2) and U(2) are irreducible tensors 
with components 

1 /Att\ ^ 
1 /47r\ ^ 

S.{2q) = -i^-) SfY,,{S,)- 2 = 1,2 (31) 

and 



t/(2g)= (^Vr2,(f). (32) 



5 J 

In dSH) and (jS2I), functions of the type 1^2 g(x) are spherical harmonics and exphcit 
expressions for |T(2)| |5) and (r||U(2)||/) are given in (US)) - (IX5|l . 



5. Spin-polarized metastable helium atoms 

For the case of metastable helium atoms, A = 0, Si = l;i = 1,2 and so the adiabatic 
potentials required for the and molecular states are the potentials Voo(r) 
and Vo2{r). Initially the atoms are spin polarized so that 5 = 2 and Ms = 2 and 
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then colhsions take place that produce final states with S' = 2,Mg = 0,±1,±2 and 
5" = 0, M'g = 0. For the state we use the analytical potential of Starck and Meyer 
[in] which has a scattering length of 157ao and supports 15 bound states. For the 
potential, we use the results obtained by Miiller et al [17 for r < 12ao and for 
r > 12ao, the potential is matched smoothly onto the long-range form Vo2(r) — Vey^chi^) 
where Vexchl?^) = Aexp(— 7r) [THl CH]- If the spin polarization is destroyed, there is 
a high probability of Penning and associative ionization and subsequent loss of atoms 
from the trap. We model this loss by using a complex optical potential of the form 
yK'S'ij) = ^00 (^) ~ |iroo('")- Two forms for roo(?") are used; rM(?") a least squares fit 
to the tabulated results in ^.17J and the simpler form rGMs('") = 0.3 exp(— r/1.086) of 
Garrison et al [20] which decreases more rapidly as r increases and does not decrease 
for small values of r. 

For this case, the change in energy of the states with / = is obtained from (jSHl), 
(EHl) and (ESI, i.e. 

AE2 = -(3^ / Ffcoo2(r)- [Gk2m{r) + 7G,202(r)] dr . (33) 

In ^21; it was shown that since the effective range of the bound-state wave function is 
typically lO'^ao to lO^ao, the wave function F^if^s for a trap state could be replaced by 
a free-wave function of the same energy and an excellent value for the energy obtained. 
Therefore in this application, the energy shifts and widths are calculated using both 
bound-state (Atrap = 1) and free- wave solutions (Atrap = 0) of (jl} for Ffcoo2 to test 
further the validity of the free-wave approximation. Energy shifts and widths are also 
calculated for trap states with 1 = 2. 



6. Numerical calculations 



The unperturbed eigenvalue equation in the form (jU or ((7j) was solved using the two 
computational methods described in ^3!- The first combines the use of quantum defect 
theory, numerical integration and an iterative procedure (QDT) and in the second a 
direct numerical solution is obtained using a discrete variable representation (DVR) 
of the kinetic energy operator and a scaled radial coordinate grid. The DVR method 
is easily modified to solve for the perturbed functions GkiAsip)- general real 
invertible transformation of the radial variable p given by 

t = u{p); p = u-\t) = U{t) (34) 

is introduced so that (pH) becomes 

d2 



Git) = Wit) Fit) 



(35) 



where 



fit)^ 



dU 



1/2 

_dt J ' 

G kv A' s'[U i^)] 

W) 



Fit) 



MAS 



Pit)] 



Wit)^ 



fit) ' 

w[uit)] 

fit) 



(36) 
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and 



Qit) 



+ 1) 
p2 




+ f 



-An* . 



(37) 



dt2 



The DVR is constructed by using a finite set of basis functions {(j)m{t)} and coordinate 
points {tm} over the interval [tiji^], so that the differential equation is transformed 
into the matrix eigenvalue equation 



where z = 1, 2, . . . A^. The matrix element Tij of the kinetic energy operator T = — d^/dt^ 
obtained using a Fourier basis is given in |13j and we choose pi^ = 2ao and pat = 15. 
The scaling is given by 



and we choose ( = 20 and p = 10 so that about 17% of the scaled mesh points lie 
between pi and ( where the interatomic and spin-dipole interactions are significant. 
Four to five digit convergence is obtained for the perturbed energies with = 2000. 

Alternatively, having determined the values of n* for the bound states, (j3)) and (fT^ 
are solved for the correct energies with Atrap = 0. The previous numerical procedure 
QDT can be readily modified for this purpose and equations (H)) and (fT^ are integrated 
numerically using the Numerov algorithm. The solution of (fT^ contains, in general, a 
particular integral plus a complementary function which is some linear combination of 
F^{r) and Fl,{r) as defined by (PT|) with IAS replaced by I'A'S'. The method, to be 
labelled as QDTF, is as follows. Equation is integrated outwards from the origin 
to some r = rmax and F^{r) and dF^{r)/dr are obtained. On matching to solutions 
that are asymptotically plane waves using methods similar to those described in |18j . 
we obtain 6i'{k). Then -F//(r) and dF;^(r)/dr are calculated at r = rmax, c.f. fj?T|l . so 
that (IH) can be integrated inwards to obtain the irregular solution. Integration of ()19p 
requires some care in dealing with the multiples of F^{r) and Fl,{r) that build up in 
the Numerov integration in both the outward and inward directions. In the integration 
of (fTTHl outwards, multiples of F^{r) are removed at each integration step determined 
so that the function at the current point is zero and this specifies G'out('")- The function 
Gin(r) is obtained by integrating inwards and a multiple of Fl,{r) is subtracted at the end 
to make the function zero at the innermost point. Finally, a multiple of Ff-{r) is added 
to G'out('^) so that the solution matches G^^ir) at r ^ 6ao and this procedure is very 
insensitive to the precise choice of matching point. This completes the specification 
of the solution Gk'i'A's'{f) of the inhomogenous equation (fT^ . It is found that for 
'"max = lOOOao convergence in the evaluation of the radial integrals in ()33|1 has been 
obtained correct to three significant figures. 



N 



J2 [f\t,)T,,f\t,) + Qit,) 5,,] Git,) = Wit,) Fit,) , 



(38) 





(39) 
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7. Results and discussion 

The scaled energy shifts and widths An* for the five lowest trap states with / = have 
been calculated using the DVR method for trapping frequencies v = ui/2tt ranging from 
1 kHz to 10 MHz. Calculations have also been carried out for frequencies 100 kHz, 1 
MHz and 10 MHz using the QDTF approximation and the agreement with the DVR 
results is very satisfactory. Differences for transitions to the S" = states range from 
3% for the lowest state to 0.4% for the highest, and at 10 MHz are always < 1%. For 
transitions to the S' = 2 state the differences are even less, they are < 0.3% except 
for the lowest state at 100 kHz where the difference is 2%. It is clear why the QDTF 
approximation is not so good for the lowest states as the range of r over which the 
bound-state wave function can be approximated closely by a free wave decreases with 
decreasing energy. The QDTF approximation has also been used in to investigate 
widths and shifts for / = 2 and they prove to be negligible. This is not surprising as 
the perturbation weights the inner region strongly where the initial wave functions with 
/ 7^ are very small. 

The results shown in table 1 have been obtained using the DVR method. The 
contributions to An* of the terms with 5" = 0, 2 in (jH!^ are also shown, together with 
the effective quantum number n* and the lifetime r given by 

T = -l/{ATTu^An*) . (40) 

The contributions from S" = are sensitive to the form used for roo('")- In table 1 we 
also show lifetimes obtained using both rM(r) and rGMsl?")- The energy shifts increase 
with trapping frequency and n*, but the fractional shifts An* /n* decrease with n* and 
are of the order 10~^ to 10^^. The energy shift arises predominantly from the transitions 
to S" = 2 states. As expected, the lifetimes tgms are shorter than the tm, as rGMs(?^) 
is larger than rM(T') for all r < Soq. This dependence suggests that an experimental 
study of lifetimes as a function of trap frequency could yield an improved knowledge of 
the potential representing the decay channel. 

The decay widths and lifetimes depend strongly on trapping frequency and the 
lifetimes range from the order of 10^ s for 1 kHz traps to the order of 15 ms for 10 
MHz traps. The lifetimes for the lower frequencies are greater than both the metastable 
helium lifetime of 8000 s and the typical lifetimes of experimental traps which are of the 
order of seconds Over the range of frequencies investigated the real and imaginary 
parts of An* can be quite closely fitted to the analytic form 

An* = -A^/{1 + Bx + Cx^ + Dx^)] x = un* , (41) 

where u is in units of MHz and 

A = 4.0144 X 10"^ 5 = 1.8806 X 10"^ 

C = 3.0740 X 10"^ ; D = -3.3204 x 10"^ (42) 
for 3?An*. For '^An* 

A = 1.7294 X 10"^; 5 = 1.7145 x 10~^ 
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C = 3.3987 X IQ-^ ; D = -3.7630 x 10"^ (43) 

and the parametric fit to the hfetimes is then directly obtained from ()4()j) . 1)411) and ()43)1 . 

In conclusion, second-order perturbation theory with a rather general form for the 
perturbation and valid for states with non-zero angular momentum has been developed 
here. These results are important and have been obtained so that the analysis can be 
directly applied to other perturbations in future. 



Appendix 

The book on Angular Momentum by Edmonds 21 is the basic reference for the notation, 
equation numbers and tables quoted below. The Wigner-Eckhart theorem, see (5.4.1), 
defines the reduced matrix elements (j'||X(A)||j) using the relation 

(/m;.|X(Ag)|jm,) = (-l)^'-^(/||X(A)||j)(_{^, ! i)' (^.1) 

where ^ ^^isa3 — j symbol and the average over the degenerate initial states 
\lm) is carried out using 

(2j + l)-i 5^ f ^' , \ ^)(^\ ^ M =[(2j + l)(2/ + l)]-iV.> (A.2) 

m'm, -' J / ^ J J ' 

J ■' 

see (3.7.8). The states \SMs) and \SiMs^) are eigenvectors of the operators S^, S^, Sf, 
and 5*,^, i = 1,2. Therefore using (jA.ljl . (4.6.3) and table 2, the reduced matrix elements 
required are given by 

(5;5^5'||5(2)||5i525) = ^ [(25 + 3)(25 + 2)(25 + 1)(25)(25 - 1)]^ 6s's ; 

{S',\S,{2)\\S,) = \ [{2S, + 3)(25, + 2){2S, + 1)(2S',)(25, - 1)]^ ; ^ = 1, 2 (A.3) 
and 

'V 2 I 


Finally, using (7.1.7) and (7.1.8), the reduced matrix elements of the spin operators Si 
and 1S2 in the coupled representation are given by 

(5;5^5'||5i(2)||5i525) = (-l)^i+^^+^[(25+l)(2y + l)]^|^i II ^2 

X {S[\\Si{2)\\Si)5s',s,; 
(5;5^5'||52(2)||5i525) = (-l)^^+^^+^'[(25 + l)(25' + l)]^|| II 2' 

X {S'2\\S2{2)\\S2)5s'^s,, (A.5) 
where |^ ^ ^|isa6 — j symbol. 

The following important special cases are easily derived from (jA.ljl - ()A.5)1 and 
tables 2 and 5 : 

|(0||T(2)||2)p = 10; |(2||T(2)||2)p = 70; | (2| |U(2)| |0) p = 1 . (A.6) 



|(r||U(2)||0P = (2/ + l)(2r + l) 



(A.4) 



Collisional loss of atoms from harmonic traps 



11 



References 

[1] Greiner M, Mandel O, Esslinger T, Hansch T W and Bloch I 2002 Nature 415 39-44 

[2] Miiller-Seydlitz T, Hartl M, Brezger B, Hansel H, Keller C, Schnetz A, Spreeuw R J C, Pfau T 

and Mlynek J 1997 Phys. Rev. Lett. 78 1038-41 
[3] Brennen G K, Caves C M, Jessen P S and Deutsch I H 1999 Phys. Rev. Lett. 82 1060-3 
[4] Jaksch D, Briegel H-J, Cirac J I, Gardiner C W and ZoUcr P 1999 Phys. Rev. Lett. 82 1975-8 
[5] Tiesinga E, Williams C J, Mies F H and Julienne P S 2000 Phys. Rev. A 61 063416 
[6] Bolda E L, Tiesinga E and Julienne P S 2002 Phys. Rev. A 66 013403 
[7] Blume D and Greene C H 2002 Phys. Rev. A 65 043613 

[8] Pereira Dos Santos F, Leonard J, Wang Junmin, Barrelet C J, Perales F, Rasel E, Unnikrishnan 

C S, Leduc M and Cohen- Tannoudji C 2001 Phys. Rev. Lett. 86 3459-62 
[9] Robert A, Sirjean O, Browaeys A, Poupard J, Nowak S, Boiron D, Westbrook C I and Aspect 
A 2001 Science 292 461-4 
[10] Sirjean O, Seidelin S, Viana Gomes J, Boiron D, Westbrook C I, Aspect A and Slilyapnikov G 

V 2002 Phys. Rev. Lett. 89 220406 
[11] Seidelin S, Sirjean O, Viana Gomes J, Boiron D, Westbrook C I and Aspect A 2002 e-print 
arXiv:cond-mat/0211112 

[12] Peach G, Whittingham I B and Beams T J 2002 Interactions of Cold Atoms and Molecules ed P 
Soldan, M T Cvitas, J M Hutson and C S Adams (Daresbury: Collaborative Computational 
Project on Molecular Quantum Dynamics (CCP6)) p 85-8 

[13] Peach G, Whittingham I B and Beams T J 2004 e-print arXiv:physics/0212003, Phys. Rev. A to 
be published 

[14] Dalgarno A and Lewis J T 1955 Proc. R. Soc. A 233 70-4 

[15] Abramowitz M and Stegun I A 1965 Handbook of Mathematical Functions (New York NY: Dover) 
[16] Stiirck J and Meyer W 1994 Chem. Phys. Lett. 225 229-32 

[17] Miiller M W, Merz A, Ruf M-W, Hotop H, Meyer W and Movre M 1991 Z. Phys. D: At., Mol. 
Clusters 21 89-112 

[18] Venturi V, Whittingham I B, Leo P J and Peach G 1999 Phys. Rev. A 60 4635-46 
[19] Leo P J, Venturi V, Whittingham I B and Babb J F 2001 Phys. Rev. A 64 042710 
[20] Garrison B J, Miller W H and Schaefer H F 1973 J. Comput. Phys. 59 3193-9 
[21] Edmonds A R 1974 Angular Momentum in Quantum Mechanics 2nd edn (Princeton NJ: 
Princeton University Press) 



Collisional loss of atoms from harmonic traps 
Tables and table captions 



12 



Table 1. Effective quantum numbers n* = E/{2hoj), shifts An* (in units of 10~^) and 
state lifetimes tm for harmonic traps with frequencies ranging from 1 kHz to 10 MHz. 
Also shown are the lifetimes tgms for comparison. Numbers in parentheses denote 

l)()\v(>rs of 10. 



< 


n* 


An*{S' 


= 0) 


An*{S' = 2) 


An* (Total) 


tm(s) 




T'GMs(s) 


i/ = 1 kHz 





0.7521 


-0.0007 - 


0.0005i 


-0.0105 


-0.0112 - 


0.0005i 


1.65(4; 




1.42(4) 


1 


1.7531 


-0.0010 - 


0.0007i 


-0.0158 


-0.0168 - 


0.0007i 


1.10(4; 




9.47(3) 


2 


2.7539 


-0.0013 - 


0.0009i 


-0.0197 


-0.0210 - 


0.0009i 


8.83(3; 




7.58(3) 


3 


3.7546 


-0.0015- 


O.OOlli 


-0.0229 


-0.0244 - 


O.OOlli 


7.57(3; 




6.50(3) 


4 


4.7551 


-0.0017- 


0.0012i 


-0.0258 


-0.0275 - 


0.0012i 


6.73(3; 




5.78(3) 


z/ = 10 kHz 





0.7566 


-0.0022 - 


0.0015i 


-0.0335 


-0.0356 - 


0.0015i 


5.20(2; 




4.47(2) 


1 


1.7599 


-0.0032 - 


0.0023i 


-0.0500 


-0.0532 - 


0.0023i 


3.48(2; 




2.99(2) 


2 


2.7624 


-0.0040 - 


0.0029i 


-0.0623 


-0.0663 - 


0.0029i 


2.79(2; 




2.39(2) 


3 


3.7644 


-0.0047- 


0.0033i 


-0.0726 


-0.0773 - 


0.0033i 


2.39(2; 




2.05(2) 


4 


4.7662 


-0.0053 - 


0.0037i 


-0.0816 


-0.0869 - 


0.0037i 


2.13(2; 




1.83(2) 


u = 100 kHz 





0.7711 


-0.0069 - 


0.0049i 


-0.1073 


-0.1142 - 


0.0049i 


i.62(i; 




1.39(1) 


1 


1.7815 


-0.0102 - 


0.0073i 


-0.1588 


-0.1691 - 


0.0073i 


i.09(i; 




9.37 


2 


2.7893 


-0.0127- 


0.0091i 


-0.1971 


-0.2098 - 


0.0091i 


8.78 




7.54 


3 


3.7957 


-0.0147- 


0.0105i 


-0.2288 


-0.2435 - 


0.0105i 


7.56 




6.50 


4 


4.8014 


-0.0165 - 


0.0118i 


-0.2562 


-0.2727 - 


0.0118i 


6.75 




5.80 


i/ = 1 MHz 





0.8181 


-0.0225 - 


0.0160i 


-0.3493 


-0.3718 - 


o.oieoi 


4.96(- 


1) 


4.26(-l) 


1 


1.8502 


-0.0318- 


0.0229i 


-0.4960 


-0.5277- 


0.0229i 


3.48(- 


1) 


2.99(-l) 


2 


2.8738 


-0.0382 - 


0.0275i 


-0.5967 


-0.6349 - 


0.0275i 


2.89(- 


1) 


2.48(-l) 


3 


3.8930 


-0.0430 - 


0.0311i 


-0.6740 


-0.7170 - 


0.0311i 


2.56(- 


1) 


2.19(-1) 


4 


4.9096 


-0.0469 - 


0.0340i 


-0.7362 


-0..7831 - 


- 0.0340i 


2.34(- 


1) 


2.01(-1) 


u = 10 MHz 





0.9703 


-0.0645 - 


0.0469i 


-1.0145 


-1.0791 - 


0.0469i 


1.70(- 


2) 


1.46(-2) 


1 


2.0509 


-0.0720 - 


0.0533i 


-1.1468 


-1.2188 - 


0.0533i 


1.49(- 


2) 


1.28(-2) 


2 


3.1040 


-0.0732 - 


0.0550i 


-1.1783 


-1.2516 - 


0.0550i 


1.45(- 


2) 


1.24(-2) 


3 


4.1442 


-0.0725 - 


0.0551i 


-1.1777 


-1.2502 - 


0.0551i 


1.44(- 


2) 


1.24(-2) 


4 


5.1767 


-0.0710 - 


0.0546i 


-1.1648 


-1.2358- 


0.0546i 


1.46(- 


2) 


1.25(-2) 



